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We analyze doping dependence of the spin resonance of an s superconductor and its sensitivity 
to the ellipticity of electron pockets, to magnetic and non-magnetic impurities, and to the angle 
dependence of the superconducting gap along electron Fermi surfaces. We show that the maximum 
intensity of the resonance shifts from commensurate to incommensurate momentum above some 
critical doping which decreases with increasing ellipticity. Angle dependence of the gap and partic- 
ularly the presence of accidental nodes lowers the overall intensity of the resonance peak and shifts 
its position towards the onset of the particle-hole continuum. Still, however, the resonance remains 
a true S— function in the clean limit. When non-magnetic or magnetic impurities are present, the 
resonance broadens and its position shifts. The shift depends on the type of impurities and on the 
ratio of intraband and interband scattering components. The ratio £l res /T c increases almost linearly 
with the strength of the interband impurity scattering, in agreement with the experimental data. 
We also compare spin response of s and s ++ superconductors. We show that there is no resonance 
for s ++ gap, even when there is a finite mismatch between electron and hole Fermi surfaces shifted 
by the antiferromagnetic momentum. 

PACS numbers: 74.20. Mn, 74.20.Rp, 74.25. Jb, 74.25.Ha 



I. INTRODUCTION 

The relation between unconventional superconductiv- 
ity and magnetism is one of the most interesting topics in 
condensed-matter physics. In conventional phonon- me- 
diated s— wave superconductors (SCs), the SC gap A is 
approximately a constant along the Fermi surface (FS), 
and paramagnetic spin excitations at T << T c are sup- 
pressed below 2A due to formation of Cooper pairs with 
the total spin S = 0. In unconventional SCs, such as 
layered cuprates or some heavy fermion materials, the 
pairing symmetry is d— wave, and the SC gap changes 
sign along the Fermi surface. In this situation, the sup- 
pression of the spin response is only one effect of super- 
conductivity, another is the appearance of the resonance 
below the onset of the particle-hole continuum, at a set of 
momenta which connect FS points with different signs of 
the gapi The intensity of the resonance is the strongest 
at Q for which the gap magnitude is the largest at the 
corresponding FS points. For many unconventional SCs 
such Q coincides with the antiferromagnetic momentum 

(7T,7T). 

In general, the spin resonance has contributions from 
both fermionic excitations in particlc-holc and particle- 
particle channels, which are mixed in a superconductor^ 
In most cases, however, the dominant contribution to 
the resonance comes from particle-hole channel and from 
this perspective the resonance can be viewed as a spin 
exciton 3 -. The spin resonance is a true S— function in a 
2D superconductor at T = 0, but acquires a finite width 
at a finite T in 2D, and even at T = in 3D systems^, if a 
line of nodes intersects the locus of FS points separated 
by Q. A broadened resonance survives even when the 



system loses superconducting phase coherence^, i.e., it 
still exists in some T range above T c , however its intensity 
sharply increases only below T c . 

Because spin resonance only develops when there is a 
sign change of the superconducting gap between kp and 
kf+Q, it is widely regarded as a probe of unconventional 
gap symmetry or structure, complementary to phase sen- 
sitive measurements. The observation of the spin reso- 
nance in the cuprates^ and in some heavy-fermion mate- 
rials— ^ is a strong evidence of a d— wave gap symmetry 
in these materials. 

The subject of this paper is the analysis of the several 
features of the spin resonance seen by inelastic neutron 
scattering (INS) in Fc-bascd superconductors (FcSCs). 
These systems are multi-orbital/multi-band metals with 
two or three hole FS pockets centered at the T point (0, 0) 
and two elliptical electron FS pockets centered at (ir, tt) in 
the actual, folded Brillouin zone (BZ) with 2 Fe— atoms 
per unit cell. There also exist FeSCs with only hole or 
only electron FSs, but we will not discuss these materi- 
als. Throughout this paper we assume that FeSCs can 
be reasonably well approximated by 2D models and ne- 
glect the anisotropy of electron dispersion along k z . We 
also consider only interactions that conserve momentum 
in the unfolded BZ (1 Fe— atoms per unit cell) and ne- 
glect additional interactions via a pnictogen, which only 
conserve momentum in the folded BZ and hybridize the 
two electron pockets. This hybridization is relevant to 
systems with only electron FSs^, but does not look to be 
important for our consideration. 

Angle Resolved Photoemission Spectroscopy (ARPES) 
studies of FcSCs with hole and electron pockets demon- 
strated quite convincingly that the SC gap in these sys- 
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terns does not have nodes along the hole FSs, at least 
for k z probed by ARPES. The absence of the nodes is 
the proof that FeSCs are s— wave superconductors, oth- 
erwise the gap would have nodes on hole pockets for all 
k z . Still, however, in a multi-band SC, s— wave gap can 
be unconventional and give rise to the spin resonance in 
INS. Indeed, s-wave symmetry only implies that the SC 
gap is approximately a constant along hole FSs, but does 
not impose a restriction on the relative signs of the SC 
gaps along hole and electron FSs separated by Q. A gap 
structure consistent with ARPES can be either a conven- 
tional, sign-preserving s ++ gap, or an unconventional s ± 
gap which changes sign between hole and electron FSs. 

Both gaps have been proposed for FeSCs based on the 
two different assumptions about the interplay between 
intra-orbital and inter-orbital screened Coulomb interac- 
tions, s gap has been proposed based on the assumption 
that intra-orbital interaction is stronger than the inter- 
orbital one. In the band description, this condition im- 
plies that both intra-pocket and inter-pocket interactions 
are positive, i.e., repulsive. Then a conventional s— wave 
superconductivity is impossible, but s ± pairing is possi- 
ble if the inter-pocket is larger than the intra-pocket one. 
A positive inter-pocket interaction is enhanced by anti- 
ferromagnetic fluctuations and close enough to a mag- 
netic instability overcomes intra-pocket repulsior>i2r— . 
The same antiferromagnetic fluctuations also enhance a 
d-wave pairing component, but a d-wave gap has nodes 
on the hole FSs and has a smaller condensation energy, 
at least near a magnetic transition. 

The alternative, s ++ pairing has been proposed^ 
based on the opposite assumption that intra-orbital in- 
teraction is weaker than the inter-orbital one. In the 
band description, this implies that intra-pocket interac- 
tion is repulsive (positive), but inter-pocket interactions 
is attractive (negative). Then a conventional s ++ super- 
conductivity occurs once inter-pocket interaction exceeds 
the intra-pocket one. A negative inter-pocket pairing in- 
teraction is enhanced by charge (orbital) fluctuations in 
combination with phonons, that is, if orbital fluctuations 
play the dominant role, the system develops an s ++ su- 
perconductivity 

Because s ± gap changes the sign between hole and elec- 
tron FSs, it satisfies the same condition for the resonance, 
Ak = — Ak+Q as in a d— wave superconductor. Accord- 
ingly, the spin response of an superconductor below 
T c should contain a sharp, nearly 5— functional spin res- 
onance^—. No such resonance develops if the gap is a 
conventional s ++ structure. 

The peak in the dynamical spin susceptibility has been 
observed below T c at the antiferromagnetic wave vector 
Q in several FeSCs— ~— and has been regarded by many 
as the strong, yet indirect experimental evidence in favor 
of gap. [Other indirect evidences for s ± SC are the ob- 
servation of a magnetic field dependence of the quadripar- 
tite interference peaks in STM2£ and the very presence of 
the co-existence phase between antifcrromagnctism and 
SC2£]. These experimental findings are particularly im- 



portant because direct phase sensitive measurements of 
the gap in FeSCs are lacking. 

Another interpretation of the magnetic peak was put 
forward by the supporters of s ++ pairing^. They argued 
that, even in an s ++ superconductor, there is a redistri- 
bution of a magnetic spectral weight below T c due to the 
opening of a spin gap. This by itself leads to the de- 
velopment of a peak in the differential dynamical spin 
susceptibility (the one below T c minus the one at T c ). 
This peak is not a resonance and occurs above 2A, con- 
trary to the resonance which is a true bound state and as 
such must be located below 2A. In principle, this differ- 
ence should be sufficient to determine which scenario is 
consistent with the data, particularly given that the gaps 
have been measured by ARPES. The difficulty, however, 
is that FeSCs have several gaps of different magnitudes, 
and the measured position of the peak in the magnetic 
response is below 2A maa: but above 2A m „ and hence can 
be interpreted both ways. 

In this paper, we analyze in detail other properties of 
the spin resonance peak in FeSCs which could help es- 
tablish more accurately whether or not the development 
of the peak in the spin response below T c is the evidence 
for s ± gap. For this, we explore in detail similarities 
and differences between the resonance in FeSCs and the 
resonance in d-wave cuprate superconductors. 

The spin resonance in the cuprates has been studied 
in great detail over the last 15 years. The FS in the 
cuprates is large and both kp and kp + Q are on the 
same FS sheet (these special points are called hot spots). 
The velocities at the two hot spots separated by Q arc 
neither parallel not antiparallcl, such that the original 
FS and the shadow one, with k shifted by Q, simply 
cross at a hot spot. In this situation, the sign change of 
the gap between the two hot spots separated by Q im- 
plies that the imaginary part of the bare particle- hole sus- 
ceptibility Imxo(Qj^) jumps at 2A from its zero value 
below 2A to a finite value immediately above 2A. By 
Kramers-Kronig relation, the real part of the particle- 
hole susceptibility Rcxo(Qi^) then diverges logarithmi- 
cally at 2A. Because Rexo(QjO) changes only little be- 
tween the normal and the SC state and Rexo(Qj^) di- 
verges upon approaching 2A from below, the full sus- 
ceptibility %(Q,fi) oc (1 — Uxo(QM))^ 1 develops a 
S— functional resonance at some O res (Q) < 2 A, where 
Imxo(Qj^) = 0. Upon momentum deviations from Q, 
Xo(cj+Q, 0) decreases chiefly due to the momentum vari- 
ation of the d-wave gap. The resonance frequency ap- 
proaches zero at q + Q equal to the distance between the 
two opposite points on the FS where the d— wave gap 
vanishes. [At even larger q the new, second resonance 
develops and the resonance frequency bounces back^.] 

In FeSCs, the physics behind the resonance is similar to 
the one in the cuprates in the sense that the sign change 
of the gap between hole and electron FSs again gives rise 
to the divergence of the particle- hole susceptibility at 2 A. 
But there are also two important differences. First, an 
s— wave gap does not induce a downward dispersion of 
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the resonance. Second, the resonance in FcSCs is more 
likely to become incommensurate either upon doping or 
even at zero doping. This happens because the geome- 
try of electron and hole FSs in FeSCs is such that there 
exist special momenta q ^ 0, at which a hole FS shifted 
by q + Q just touches an electron FS. At these special q, 
the real part of the particle-hole susceptibility diverges as 
(2A — SI) -1 / 4 , i.e. stronger than logarithmically. These 
q appear at a finite doping if all FSs are approximated 
by circles and exist already at zero doping if we treat 
electron FSs as elliptical. Because the stronger is the di- 
vergence, the farther down is the resonance from 2A, the 
position of the lowest f2 res shifts, upon increasing doping 
or ellipticity, from a commensurate Q to an incommensu- 
rate q+Q, with a non-zero q. We show that the intensity 
of the resonance does not follow this shift instantly, but 
eventually the intensity also becomes the strongest at an 
incommensurate q+Q. A transformation from a com- 
mensurate to an incommensurate spin resonance upon 
doping has been found in the recent neutron study of 
the evolution of the resonance in Bai_ 2 -K a; Fe2As2 with 
doping^. An incommensurate spin resonance has been 
observed in Refi^ in a superconducting FeTco.6Seo.4i2 1 . 

We also discuss in detail the width of the resonance 
peak in FeSCs. Experimentally, the peak is quite broad 
already at low T, which was argued^ to be inconsis- 
tent with a true resonance. We show that there are two 
effects, specific to FeSCs, which lead to additional broad- 
ening of the peak. First is the symmetry-imposed cos 29 
angle variation of the gap along electron FSs, the sec- 
ond is pair-breaking effect on superconductivity from 
non-magnetic impurities which scatter between hole and 
electron FSs. 

We show that the impurity scattering broadens the 
resonance peak, but its position, i.e. f2 res , remains es- 
sentially unaffected. Given that T c decreases by impurity 
scattering, the ratio £l res /T c increases with increasing im- 
purity concentration which correlates with the increase 
of doping. This trend is in agreement with recent INS 
experiments^ which show that Q res /T c moves up with 
increasing doping. 

The two key conclusions of our study is that (i) the 
doping dependence and the shape of the spin resonance 
in FeSCs is fully consistent with the gap structure, 
and (ii) the spin response of an SC can be fully un- 
derstood within an itinerant approach, without invoking 
any localized moments. 

The paper is organized as follows. In Sec. [II] we intro- 
duce the model and discuss the computation of a spin 
susceptibility in a multi-band superconductor. In Sec. 
IIIII we discuss the doping dependence of the resonance, 
doping-induced incommensuration, the role of ellipticity 
of the electron pockets, and the role of the angle varia- 
tion of the gap. In Sec. HVI wc discuss how the resonance 
is affected by non-magnetic impurities, which are pair- 
breaking for s gap. Results with magnetic impurities 



are also included. Our conclusions are presented in Sec. 

For completeness, in the Appendix we discuss one 
particular aspect of the interplay between magnetic re- 
sponses at Q in s ++ and superconductors. Namely, 
previous calculations have demonstrated that there is no 
resonance in an s ++ superconductor, when the original 
and the shadow FSs cross. We verify what happens in 
an ,s ++ SC if the original and the shadow FSs do not 
cross, which is the case when both hole and electron FSs 
are near-circular and the doping is finite, such that one 
FS is larger than the other. In this situation, the onset 
of particle hole continuum in Imxo(Q, shifts upwards 
from 2A. We find that, even in this case, there is no 
resonance for s ++ gap (and there is the resonance for s 
gap). The only feature that emerges for s ++ gap is a 
broad maximum in Jmx(Q,r2) above 2A. 



II. THE MODEL 

We start from the 4-band model with two circular hole 
pockets at (0, 0) (a-bands) and two elliptic electron pock- 
ets at (tt, 0) and (0, tt) points in the unfolded Brillouin 
Zone (UBZ) (/3-bands). The quadratic part of the Hamil- 
tonian is 

p,er,i— 1,2 

To facilitate numerical calculations, we consider lat- 
tice dispersions for all four bands, although the reso- 
nance essentially comes only from fermionic states near 
the FSs. We set e"* = t a (cosp x + cosp. y ) — [ii and 
e^ 1 = e + tp([l + e] cos(p x + 7r) + [1 - e] cos^)) - Mi, 

£ p = e o + tp (t 1 - £ ] cos(p x ) + [1 + e] cos(p y + tt)) - /xl 
The momenta are measured in units of inverse Fe-Fe 
spacing a x = a y = a. The parameter e accounts for 
the ellipticity of the electron pockets. To make qualita- 
tive comparisons to experiments, we set Fermi velocities 
and Fermi surfaces to be equal to those in Refsi 32 ' 33 For 
this we take t a = OMeV, tp = -0.68eV, m = 1.54eV, 
= 1.64eV, and e = 0.31eV. For e = 0.8 the Fermi 
velocities are 0.5eT^a for the two degenerate a-bands, 
and v x = 0.27eVa and v y = Q.AQeVa along x- and y- 
directions for the /3i-band (and vice versa for the fa 
band). The Fermi surfaces are shown in Fig. Ufa). 

The interacting part of the Hamiltonian contains four- 
fcrmion interactions with small momentum transfer and 
momentum transfers (tt, 0), (0, tt), and (tt,tt). They in- 
clude a — (3 interactions with momentum transfer (0, tt) 
and (tt, 0) as well as /3-/3 and a-a interactions with mo- 
mentum transfers (0,0) and (tt, tt), respectively, and are 
given by 
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FIG. 1: color online) (a) Calculated Fermi surfaces for the 
four band model. Solid and dashed lines for the electron FSs 
are for e = 0.5 and e = 0.8, respectively. The arrow shows 
the deviation from perfect nesting due to ellipticity; (b) Dia- 
grammatic representations of the interactions in the four band 
model. 



The vertices are shown diagrammatically in Fig. [TJb). 
For simplicity, we approximate all interactions as angle- 
independent, i.e., neglect the angle dependence intro- 
duced by dressing the interactions by coherence factors 
associated with the hybridization of Fc d— orbitals. These 
coherent factors do play a role for the structure of the 
SC gap^s, but do not substantially modify the spin reso- 
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In the magnetically-disordered state transverse and 
longitudinal components of the spin susceptibility are 
undistinguishable, and we focus below on the transverse 
part. In the matrix RPA approximation, which we adopt, 
the transverse components of the full spin susceptibility 
X l ' J , where i and j are band indices, are related to trans- 
verse components of the bare susceptibility x<f as 



Xd 



,j'X 



(3) 



The summation over repeated band indices is assumed 
and u l 3 are matrix elements of the interactions shown 
in FigOJb) (contributions from all u k arc included in 
The solution of Eq. ^ in matrix form is straightforward: 
X = Xo(l - "Xo) -1 - 



The components of the bare spin susceptibility xo = 
Xo"'( < li i^m) ar e given by usual combinations of normal 
and anomalous Green's functions 

X ij (p, iO m ) = ~ J2 Tr [G*(k + p, iw„ + ifi m )G»'(k, ioj n ) 



F 4 (k + p, iu)„ + iQ m )F j (k, iw n )] 



and F l (k, ioj n ) 



where G l (k, i(J n ) = — - 

^+(^) 2 +( A L) 2 ' 

The main contribution to the full spin susceptibil- 
ity with momenta near Q comes from susceptibilities 
and interactions which involve hole and electron states 
(particle-hole bubbles made of one electron and one hole 
propagator, and u% and U3 interaction terms in Eq. ([2]). 
see RefJ£). For completeness, in numerical calculations 
we will keep all terms in the matrix equation for the full 
susceptibility. The results for the full £^ . x lJ (Pi i^m) ob- 
tained this way do not differ much from those obtained 
by keeping only electron-hole bubbles and u\ and W3 in- 
teraction terms. 

We do not present explicitly interactions leading to su- 
perconductivity and not discuss the solution of the pair- 
ing problem. This has been done in numerous other 
works on this subject^. We take the results of these 
studies as input and set the gap to be of s ± type, with 
different sign between hole and electron pockets. More 
specifically, we set the gap to be a constant along the 
hole FSs, A" 1 - 2 = A^, and set the gap along the two 
electron FSs to be A^ 1 ' 2 = A e (l ± rcos2<^>), where <fi is 
the angle counted from k x direction in the UBZ. The 
gap A' 31 * 2 has no nodes if r < 1 and has accidental 
nodes when r > 1 at non-symmetry selected directions 
cos2(/> = 1/r. The numerical results presented below are 
obtained for u\ = 1*3 rj 0.25eV ( these numbers guarantee 
that system remains paramagnetic in the normal state), 
(7) 

and 1*5 = U2 = 0.5«i, u\ = 0.1«;. For simplicity we set 
A h = -A e = A although A h and -A' 9l ' 2 (0 = tt/4) do 
not have to be equal. For the gap we used A = O.Q2t a . 
For better convergence of numerical series we added a 
small damping T = lmeV to fermionic dispersion. 
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(a) (b) (c) 

FIG. 2: (color online) Schematic form of the momentum dis- 
persion of the spin resonance frequency Cl res (solid line) at 
various doping and at zero ellipticity. Red dashed region is of 
the particle-hole continuum. The formulas show how the bare 
spin susceptibility diverges at the bottom of the continuum. 
Panel a - zero doping, panels (b),(c) - a finite doping (larger 
for c). 



III. THE RESULTS 

We present results systematically in three installations. 
We first consider the doping evolution of the resonance 
for a simple plus-minus gap and circular FSs, then we 
include the ellipticity of electron pockets, and finally we 
also include the angular dependence of the SC gap. For 
consistency with the experiments, we show all results in 
the folded BZ, when the commensurate resonance is at 
Q = (7r,7r). 

A. Doping dependence of the resonance for 
plus-minus gap and circular FSs 

This case is captured by setting e = and r = in 
A* 31 - 2 . The doping dependence is parameterized by the 
change in the chemical potential <5/i. 

Much in this case can be understood analytically. At 
zero doping, xo(Q> CI) has a strong, square-root singular- 
ity at CI = 2A. Imxo(Q,Cl) diverges as l/\/Cl — 2A at 
CI > 2A and i?ex (Q, CI) diverges as 1/y/Cl — 2A at Cl < 
2A). For incommensurate momenta, Imxoiq+Q, Cl) un- 
dergoes a finite jump at Cl = 2A and i?exo(q + Q, 2A) 
diverges logarithmically, as long as q < 2k p . In this 
situation, the lowest Cl res and the largest intensity are 
at the commensurate momentum Q. We illustrate this 
in Fig. [2£a). At a finite doping, Xo(Q>^) still have a 
square-root singularity, but now the lower boundary of 
the particle-hole continuum shifts up and the singularity 
is located at a larger CIq = \J (2A) 2 + Sfj, 2 . The reso- 
nance in the full x(Q>^) is located below CIq, but be- 
cause CIq increases with Sfi, the position of the resonance 
at the commensurate Q also increases. Meanwhile, there 
exists the range of incommensurate momenta q + Q for 
which the bottom of the particle-hole continuum is still 
located at 2Ao. These are q m in < q < qmax-, where 
q m in = \/ ' kp + <5/j 2 - k F and q max = y/ kp + Sfi 2 + k F . 
For q inside this range, 7mxo(q+Qj CI) undergoes a finite 
jump at CI — 2A and i?exo(q + Q, 2A) diverges logarith- 
mically, as before, but at the end points, i.e., at q = q m i„ 
and q = q m ax, hole and electron FSs touch each other 



after a shift by Q, and Xo(q + Q, CI) diverges by a power- 
law, this time as 1/(2A — CI) 1 / 4 . Because of stronger 
divergence, the resonance in the full x(q + Qi at these 
q shifts down from 2A more than for other q, hence the 
dispersion of the resonance develops minima at q - — Qmin 
and q = q max . 

At small doping, these two minima at incommensurate 
q are local minima because f2 res (q + Q) had an upward 
dispersion around Q at zero doping, one needs some finite 
doping to change the sign of the slope at q = 0. As a 
result, at small but finite doping, the dispersion of the 
resonance has the global minimum at q = and two 
roton-likc minima at q m i n and q m ax (see FigJ5Jb)). Once 
the doping gets larger, the resonance energy at q = 
keeps going up together with the bottom of the particle- 
hole continuum, and eventually the global minimum of 
the resonance dispersion discontinuously shifts to a finite 
q, i.e., to an incommensurate momentum (see Figj^c)). 

We reproduced this behavior in the numerical analysis. 
We display the numerical results in Fig. [3] The upper 
panels show the dispersion of the resonance energy, the 
lower panels show the imaginary part of Xo(q + Qj^)- 
FigEI a ) shows the behavior at zero doping. The 'square- 
root' divergence of imxo(Q>^) at = 2A is clearly 
visible, while for nonzero q < q ma x = 2fc p we observe 
a 'jump' in 7mxo(q + Q>^) at CI = 2 A, in agreement 
with analytical consideration. We also see, at frequen- 
cies above 2 A, the "upward" dispersion of the maximum 
in 7mxo(q + Qj CI). This dispersion originates from the 
momentum dependence of the static Xo(q+Q) and gives 
rise to the upward dispersion of f2 res (q + Q) in the upper 
panel of Fig. 03a). 

In Figs. [UJb) we show the behavior of Imxo and the 
dispersion of the resonance at small but finite dopings 
x = 0.01 (8fi = —0.02). We again see the square root 
singularity in 7mxo(q + Q, CI) at the commensurate mo- 
mentum q = 0, but now it shifts to a frequency above 2A, 
consistent with analytical CI = y/ (2A) 2 + S/j, 2 . We also 
clearly sec additional singularity in imxo(q + Q,0) at 
an incommensurate momentum, which we associate with 
Qmin (other momentum, g max is outside the momentum 
range in [3] and 7roxo(q + Q, CI) at such q already quite 
small). The resonance frequency f2 res (q + Q) still has 
a minimum at q = 0, but because of an additional sin- 
gularity in Jmxo(q + Q, CI) at q = q TO , n , the dispersion 
of Cl res around the minimum becomes rather flat. As 
doping is further increased to x = 0.03 (S/j, = —0.051) 
(FigE£c)), the square root singularity at q = is pushed 
further outward and r2 r es(Q) eventually moves to above 
2A, while the one at q = q m i n stays below 2A. As a re- 
sult, the minimum of Cl res shifts from q = to q = q m i ra , 
as is clearly seen in the upper panel of Figj3jc) . 

We found that the maximal intensity of the resonance 
does not immediately follow the position of the minimum 
of fires and over some range of dopings is still the largest 
at the commensurate momentum Q even when the mini- 
mum of Cl res already shifts to an incommensurate q+ Q. 
But as doping increases even further, the maximum of 
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FIG. 3: (color online) The imaginary part of the total physical 
susceptibility Imx(q + Q, O) for various doping levels at zero 
ellipticity: (a) x = 0, (b)x = 0.01(5/Lt = -0.02), and (c)x = 
0.03(5 /i = —0.051). The lower panel shows the corresponding 
behavior of the imaginary part of the bare spin susceptibility. 
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FIG. 4: (color online) The imaginary part of the total physical 
susceptibility Imx(q + Q, f2) for various degrees of ellipticity 
at zero doping: (a) e = 0.1, (b)e = 0.3, and (c)e = 0.5. The 
lower panel shows the behavior of the imaginary part of the 
bare spin susceptibility. 



the intensity eventually shifts to an incommensurate mo- 
mentum. 



B. Role of ellipticity of the electron pockets. 

We next consider the case of elliptical electron pockets, 
when e is nonzero. In this case, resonance may become 
incommensurate already at zero doping. Indeed, at x = 
and non-zero ellipticity, Im xo(Q> ^) has a finite jump at 
2A because hole and electron velocities at each of the 
four hot spots (points where a hole pocket shifted by Q 
crosses an electron pocket) are neither parallel nor an- 
tiparallcl. In other words, a square-root singularity at 
Q is cut by e and is replaced by a jump. However, the 
power- law, 1/(2 A — Q) 1 / 4 , singularities at incommensu- 
rate q min and q max survive. Like before, q mi „ and q max 
are the points where the hole FS, shifted by Q, touches 

nin four 

nil r • 

the 



an elliptical electron FS. At the smaller q = c 
hot spots transform into two, at a larger q = c 
remaining two hot spots disappear. Because the shift of 
fires from 2A is larger, the stronger is the divergence 
of Im xo(Qj fi) at n = 2 A, the minimum of the disper- 
sion of fl res should eventually shift to an incommensurate 
momentum once ellipticity becomes large enough to over- 
come a positive curvature of the dispersion near Q which 
persists up to some finite e. We found the same trends 
in the numerical calculations. In particular, in FigUwe 
show the evolution of the dispersion of the resonance and 
the form of Imxo(q + Q,^) a s a function of ellipticity 
at zero doping. We clearly see that once e increases, the 
dispersion becomes more flat near the minimum at Q, 
and eventually the minimum of the dispersion shifts to a 
non-zero q, i.e., to an incommensurate momentum. 

In Fig. [5]we combine finite doping and finite ellipticity. 
We see the same tends as before, namely at a given non- 
zero e (we set e = 0.3 for definiteness), the resonance is 
quite broad already at x = (Fig. [SJa)), it gets even 




X) 0.9 ! 



FIG. 5: (color online) The imaginary part of the total phys- 
ical susceptibility Imx(q+Q,fi) for various doping levels 
[(a) x = 0(5^ = 0), (b)x = 0.02(5^ = -0.035), and 
(c)x = 0.05(5/i = -0.076) and e = 0.3. The lower panel 
shows the corresponding behavior of the imaginary part of 
the bare spin susceptibility. 



more broad as x increases (Fig. EJ^b)), and eventually 
incommensurate minima appear at large enough x (Fig. 
03c)). Note that at non-zero ellipticity the development 
of incommensurate minima in £l res at is less pronounced 
effect than at e = and the key effect at small doping 
is the broadening of the resonance peak. This is because 
the bottom of the particle-hole continuum at Q does not 
move up with doping as long as an electron ellipses and a 
hole circle shifted by Q cross, and fl res at Q and at q+Q 
with q = q m in remain near-equal. Still, at larger dopings 
the resonance definitely becomes incommensurate, i.e., 
the minimum of fi res and the largest intensity shift to 
a incommensurate position q = q TO j n . There is another 
local minimum at q = q„ la2 -, like at e = 0, but that one 
corresponds to large enough q where Im\Q is already 
small. 
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C. Role of the gap anisotropy 

Another interesting aspect of the physics of the spin 
resonance in FcSCs is the role played by the angu- 
lar dependence of the s H -wave gap. We remind that 
in FeSCs the gaps on the two electron FSs behave as 
A* 912 = — A(l±r cos 20), and r is generally a finite num- 
ber. On general grounds we expect that the resonance 
gets broader simply because 2A becomes "soft" variable 
once the magnitude of 2 A varies along the electron FS. 
Only at T = and in the idealized case of no impurity 
scattering this "softness" does not matter because only 
fcrmions in the immediate vicinity of hot spots contribute 
to singularity in 7roxo(q+Q, an d these fermions have 
some fixed A& and Afc +q+ Q. At a finite T and/or in 
the presence of some residual scattering, the momentum 
range for fcrmions contributing to the resonance increases 
and the angular dependence of the electron gap becomes 
relevant. Numerical calculations arc in line with this rea- 
soning. In Fig JB] (lower panel) we show the evolution with 
r of the bare susceptibility at Q taken between a\ and 
/3i bands. Other interband susceptibilities show qualita- 
tively similar behavior. We set ellipticity and doping to 
be non-zero (e = 0.5, x = 0.03), so there is no power- 
law singularity at Q. For r = 0, Im Xo 1 ^ 1 (Q>^) then 
has a finite jump at 2A. We see that the jump persists 
also at a finite r, but its magnitude is reduced and the 
jump essentially disappears at r > 1. The reduction of 
the jump affects the behavior of the full spin suscepti- 
bility, which we show in Fig. |BJ (upper panel). We see 
that the resonance peak, which is already quite broad at 
r = 0, become even more broad at a finite r, and at large 
enough r the width of the resonance peak is essentially 
controlled by r. 

The broadening of the peak at Q is wider than in the 
cuprates because there the FS is large and the momenta 
which contribute to the resonance arc confined to near 
vicinity of hot spots, where Afc and A^+q arc rather flat 
and can be treated as constants^. In FeSCs, the FSs are 
much smaller and, accordingly, the effect of gap variation 
on the width of the resonance is stronger. 

Strong enhancement of the width of the resonance peak 
due to angle variation of the SC gap is the possible ex- 
planation why the observed peak at Q is wider than one 
should expect from the fully gapped s ± -state^. Another 
possible explanation is the effect of impurities, which we 
consider in the next section. 

The effect of the gap variation on the dispersion of the 
resonance is strongest at Q because for incommensurate 
momenta and, particularly, for q = q m i n , the gaps on hole 
and electron FSs differ quite substantially by magnitude, 
and, when r > 1, have the same sign. A simple reason- 
ing shows that, in this situation, two effects happen: the 
bottom of the particle-hole continuum goes down near 
Qmin, and, at the same time, the intensity of the reso- 
nance peak at q m i ra drops, and for r > 1 the resonance at 



q ~ Qmin simply does not exist. The outcome is that, as 
r increases, the resonance dispersion flattens up near Q 




FIG. 6: (color online) The imaginary part of the bare inter- 
band susceptibility lmx ai131 (Q, w) (lower panel) and the total 
physical susceptibility (upper panel) for (a)r = 0.5, (b)r = 1, 
and (c)r = 2. We set x = 0.03 and e = 0.5. 

and the intensity drops away from Q. As shown in FigjBJ 
the numerical analysis is again in line with this reasoning. 

IV. EFFECT OF IMPURITIES 

In this section we study how the position and the width 
of spin resonance are affected by elastic impurity scatter- 
ing. We consider both magnetic and non-magnetic impu- 
rities. To keep the calculations tractable we make several 
simplifying assumptions: (i) we only consider resonance 
at T = and at the commensurate momentum Q, (ii) we 
neglect ellipticity of electron pockets and the angle varia- 
tion of the SC gap, and (iii) we approximate the band dis- 
persion by the two-band model of one hole and one elec- 
tron pockets (a- and /3-fermions) separated by Q = (tt, tt) 
and perfectly nested, i.e., set e Q (k) = — e /3 (k + Q) and, 
in the absence of impurities, A Q = — A^ = A. We show 
that impurity scattering broadens the resonance but its 
energy position f2 res (Q) can remain essentially intact. 
Because T c for pairing drops with impurity scatter- 
ing, the ratio Q res /T c increases, in agreement with the 
experimental data. Because both the broadening of the 
resonance and the increase of £l res /T c with impurity scat- 
tering are related to gap symmetry, we expect that 
these two results survive for more realistic 4-5 band mod- 
els with elliptical electron FSs and angle- dependent gap. 

A. Method 

Non-magnetic impurities affect electrons via the impu- 
rity scattering potential in the form 
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H imp = ^2 [Uq(Mi - k 2 



ki,2 



/3t ia 5**'Pk 2 v) + t^(ki-k 2 ) 



5 k 2 cr' 



(5) 



r 



where £/o(ki — It2) and J7,r(ki — k2) are intra- and inter- 
band scattering terms, respectively. For Uq, ki and k2 
are both near (0,0) or Q, while for U v ki w (0,0) and 
k 2 ~ Q or vice versa. Quite generally one expects U„ < 



J 



C/o for any impurity potential with a finite range in real 
space. 

Paramagnetic impurities in turn affect electrons via the 
spin-dependent potential 



ki. 



Hprnag = ^ Uq (kl - k 2 ) [OL^^aa' ' Sa^p + /j£ 1<T dW' ' SPkrf 

£/£(ki - k 2 ) (al ia aaa> ■ Sfik 2(7 ' + /?k lC .<W ■ Satk 



(6) 



where a = (er^, <x y , ov) and <7j are Pauli matrices, 
and S is the static impurity spin with the property 



<C SjSj ^impurity sites' 



gOg+i) 



Fermionic propagators for an s superconductors in 
the presence of non-magnetic impurities have been con- 
sidered before (see e.g., Ref. llOh . The extension to 
include paramagnetic impurities is straightforward, and 
yields 



Z iu 



^« + |A u p) + ( e ^)2 



F° 



with 



± 

Z = 

2 

b = 



ZK 



(7) 



u 



A 2 , 



w 2 „ + I A, 
2 (u^ + tig) 



(8) 



where we have adopted, for brevity, the notation u = 
ni m p\Uo(0)\ 2 (and a similar for u w , Uq, and itJJ), with ni mp 
being the impurity concentration. In Eq. ((5J, A w is the 
actual, frequency-dependent gap in the presence of im- 
purity scattering, and A& is the order parameter related 
to A u via 



As 
A 



dbJ r 



Jo ""* X 



An 



(9) 



For any 6 ^ 0, a non-zero fermionic DOS iV(a;) = 

(m/27r)i?e uj/^/uj 2 — A 2 extends to frequencies below 

2 A. For < 6 < 1, the system still preserves the gap 
in the fermionic DOS near zero frequency, i.e., over some 
range of u near cj = 0, N(lj) = 0. For b > 1, the sys- 
tem enters into a gapless regime, in which A w scales as 
at the lowest frequencies and N(lu = 0) is non-zero, 
although reduced compared to the DOS in the normal 
state. 



Observe that the gap renormalization comes from in- 
tcrband non-magnetic impurity scattering and intraband 
magnetic impurity scattering which arc pair-breaking for 
superconductor (both also reduce T c ). At the same 
time intraband non-magnetic impurity scattering and in- 
tcrband magnetic impurity scattering do not contribute 
to b and therefore are not pair-breaking. If only these two 
scattering were present, b would be zero, A w would be 
equal to A, and Eq. © would then yield A w = Aj, = A, 
i.e., the gap would not be affected by impurities. The 
same analysis extended to a finite T shows that uq and 
ut also do not affect T r . 



The next step is to use normal and anomalous Green's 
functions from ([7]) and compute spin susceptibility for 
a dirty BCS superconductor. This has to be done with 
care, though, as one has to include not only self-energy 
corrections, incorporated in (J7J, but also series of ver- 
tex corrections which are of the same order as self- 
energy terms. We follow the treatment by Gorkov and 
RusinoA*^, sum up ladder series of vertex correction dia- 
grams and after straightforward algebra obtain 
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xo(Q,ty = N J ^ 



1 + (it + u v - + u^))F A 2 + (u - + |(ufj - u^))^f 



(10) 




(c) 



b=0 



C! . 

6 

E 



FIG. 7: Calculated Imx(Q, fi) for scattering by non-magnetic 
impurities, for different values of b, which scales with the 
strength of the interband impurity scattering. We used 
uo = 1.2A (a) and uq = 1.5A (b). Figure (c) shows, for 
comparison, the resonance peak for 6 = 0. 



Resonance Peak Position 




b = 2u„/A 

FIG. 8: Scattering by non-magnetic impurities. Evolution of 
T c and the resonance peak position (both in units of A) as a 
function of b. 



where No is the 2D density of states, and 

fufu+n (fu + fu+n + 2t*o + 2u n ) 



F A 2 = 



F^a = 



wA u+S j - (w + f2)At, 



fufu+n (fu> + fu+n + 2u + 2u n ) 



(11) 



(12) 



With /„ = y/Al - ufl. 

As before, we assume that the full susceptibil- 
ity x(Q,^) is related to xo(Qj^) as x(Qj^) = 
Xo(Q,^)/(l - u apin Xo(Q,ty)- Where u spin is the resid- 
ual interaction in the spin channel (see Se dll)) . 

We solved for A w , substituted the result into (|10[) . in- 
tegrated over frequency and obtained bare Xo(Q, ^) and 
the full Im x(Qi ^) as functions of impurity scattering. 



B. Results 

We first discuss the role of non-magnetic impurities uq 
and u v and for this purpose set Uq and to zero in Eq. 
flSJ). Both uq and u n rcnormalize quasiparticle Z, but 
only contributes to pair-breaking parameter b. 

In Fig[7ta)-(b) we show the calculated Imx(^) for var- 
ious values of b and for two different uq. Imx(^) for 
b = is shown in panel (c) for comparison. We see 
that for b ~ 0.3 — 0.6 used earlier to fit NMR, pene- 
tration depth and uniform spin susceptibility data (see 
Refs i 10 i 36 ) the resonance peak in Imx(Q, f2) falls into the 
frequency range where the spin response in a dirty s^ 1 SC 
is no longer a ^-function, and its width increases with 
increasing b. At the same time, the position of the res- 
onance remains almost intact. Despite that Af, and T c 
obviously decrease when b increases. To emphasize this 
point we plot in FigjH] fl re s as a function of b together 
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FIG. 9: The effect of (a) intra-band and (b) inter-band 
impurity scattering on the resonance peak in iron-based su- 
perconductors. In particular, (a) shows the influence of the 
increasing intra-band non-magnetic and magnetic scattering 
potential, while (b) displays the corresponding effect of the 
inter-band impurity scattering. In both cases, the solid and 
dashed curves refer to the non-magnetic and magnetic im- 
purity scattering, respectively. While studying the effect of 
magnetic impurities, we kept the non-magnetic impurities as 
constants and vice versa. 



with T c (b) calculated from the conventional relatio: 



,3(5 



' 2 + 2?rT r 



(13) 



where ip is the digamma-function and T c q is the transi- 
tion temperature in the clean limit. We clearly see that 
T c decreases with b while Q res remains practically un- 
changed, such that the experimentally measured ratio 
£l res /T c increases with the strength of intra-band impu- 
rity scattering. This trend is consistent with the observed 
doping dependence of £l res /T c in Bai_ a: K. E Fc2As2 above 
optimal doping, where the decrease of T c with increasing 
K concentration is believed to be at least partly due to 
increased impurity scattering^. 

We next keep Uq and non zero and compare the 
effects of non-magnetic and magnetic impurities. 
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In Fig|9|a) we show how the resonance profile evolve 
once we increase small momentum impurity-scattering 
terms uq and Uq. We see that in both cases the peak 
broadens but it is shifted in opposite directions - up- 
wards when uq increases and downwards when Uq in- 
creases. Also, the broadening is much stronger under 
the change of Uq. FiglUJb) shows the effect on the peak 
of magnetic and non-magnetic inter-band scattering u v 
and uP. We see that the increase of u v strongly broadens 
the peak and weakly shifts it downwards, while the in- 
crease of u p tt has weaker effect on the width of the peak 
and shifts its position slightly upwards. This behavior is 
expected because, like we said, in an s ± -wave supercon- 
ductor the interband non-magnetic impurity scattering 
u n behaves much like a magnetic intra-band scattering, 
Uq. Therefore u n and itg should give rise to qualitatively 
similar effects on the resonance peak. This is indeed what 
we have found (compare arrow (2) in Fig. IHJa) and ar- 
row (1) in Figfjjb)). Likewise, the interband magnetic 
scattering, u p and the intraband non-magnetic scatter- 
ing Mo should have similar effect on the resonance, and 
this also agrees with Fig. [5] (compare arrow (1) in Fig. 
[UJa) and arrow (2) in Fig[H{b)). One has to bear in mind, 
however, that the effects of u T and u p and of uq and u% 
are similar but not equivalent because vertex corrections 
behave differently for magnetic and non-magnetic impu- 
rity scatterings (see Eq. lfTU)) '). independent on whether 
this is intra- or inter-band scattering. Still, the results 
in Fig. IHl clearly show that pair-breaking u„ and Uq have 
the strongest effect on the resonance peak. Both impu- 
rity scatterings broaden the peak but only slightly affect 
its position such that £l res /T c increases when the con- 
centration of either non-magnetic or magnetic impurities 
increases. This is the key result of this section. 



may actually be even stronger than in d— wave cupratc 
superconductors. 

We further analyzed the sensitivity of the resonance to 
non-magnetic as well as magnetic impurities. We obtain 
that both non-magnetic and magnetic impurities have 
pair-breaking components which broaden the resonance 
but do not shift much its position compared to the clean 
case. The reduction of T c by the same impurity scattering 
is much stronger, such that the ratio f2 res /T c increases 
with the strength of impurity scattering. 

The transformation of the minimum in the dispersion 
of the resonance and the maximum of its intensity from 
the commensurate to the incommensurate momentum 
upon doping is consistent with recent study of the dop- 
ing evolution of the resonance in Bai_ a: K a ;Fe2As2^. The 
increase of fl res /T c with the increase of impurity scatter- 
ing is consistent with the observed increase of this ratio 
with doping in overdoped systems, where the decrease of 
T c with doping can at least partly be attributed to the 
effect of impurities. 
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VI. APPENDIX 
VII. 8 ± VERSUS S++-WAVE GAP 



V. CONCLUSION 

In this work we have shown that the intensity and po- 
sition of the resonance peak depends on several factors: 
doping, cllipticity of the electron pockets; gap anisotropy 
and the presence accidental nodes on the electron pock- 
ets, and magnetic and non-magnetic impurity scattering. 

We find that with no cllipticity, the resonance starts 
at the commensurate momentum below 2A and moves 
outward in energy with doping. At some doping, the 
minimum in the dispersion of the resonance jumps to an 
incommensurate momentum, and at even larger doping 
the intensity of the resonance becomes the largest at the 
incommensurate momentum. With finite cllipticity, the 
intensity becomes the largest at the incommensurate mo- 
menta at smaller dopings, and for large enough ellipticity 
the resonance becomes incommensurate already at zero 
doping. Besides, at finite ellipticity, the resonance peak 
gets quite broad The inclusion of the angular variation 
of the s ± — leads to further broadening of the resonance 
peak, particularly when the gap has accidental nodes on 
the electron pockets. The broadening of the resonance 



In this Appendix, we briefly address the issue of the 
interplay between spin response in superconductors with 
and s ++ gaps in a situation when the original and the 
shadow FSs do not cross, which is the case when both 
hole and electron FSs are near-circular and the doping 
is finite, such that one FS is larger than the other. The 
case when the original FSs cross is quite similar to the 
cupratcs, and the absence of spin resonance for s ++ gap 
in this situation has been discussed earlier—. 

We show that there is no resonant enhancement of the 
Im% for the s ++ -wave case even when the original and 
the shadow FSs do not cross. To do so we look into the 
case with no impurity scattering and no cllipticity and 
compute 

X f(Q,ity = -N tJ2 I ' deG a {iu Jni e)G f \iu Jn +n,e+5 l i) 

n J 

, (14) 

where 8 a is proportional to the mismatch in FS radii 
due to doping. And a and (3 refer to the hole and elec- 
tron FSs respectively. In the absence of impurities we 
return to the case where A w is frequency independent 
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FIG. 10: (color online) Calculated frequency dependence of 
the imaginary part of the total physical susceptibility for 
(a)e = 0.3 and (b)e = for the s ++ -wave superconductor. 
We set x = 0.01 (Sfi = —0.02). We clearly see that there 
is no resonant enhancement of the full spin susceptibility for 
s ++ gap. To soften S— function in the a case we averaged 
Im x(Qi CI) around Q with Sq = ±0.1Q 



and simply equal to A. In the normal state we get 



This implies that Rexo has a square-root singularity 
at CI = 2A and hence the full Imx s (Q,f2) has a pole 
(resonance) at some CI < 2A. 

For s ++ case A = Aq, and we get 



Imx 



s++ 




N - 



2A 

7T 



(19) 



As expected, there is no jump in 7mXo ++ (Q,^) a ^ ^A 
implying no singularity in Re\o and hence no resonance 
in Im x s++ (Q, CI) below CI = 2 A. 

Consider next the case when <5/z is non-zero. Now 
lmxo{Q,Ci) has a singularity at CI — \J S/j, 2 + 4A 2 . Ex- 
panding around this CI by introducing a small z such that 

(§) 2 -(^ ^ 



A 2 + z 2 , we find using Eq. [TH that 



Im Xo (Q,Cl)=N —6(Cl-\5 f i\) 



(15) 



In the SC state (after appropriately adding the F*F 
term), we get 

7T 2 ^ A(A + A Q ) + (w + x)(Cl - 5n) 



Im Xo (Q,n) = -N —Y, 



\Clx — S/j, to\ 



(16) 

where A and Aq are the gaps on the hole and electron 
FSs, and the two values of x which are summed over are 
given by 



5fi CI 
Y 2 



(t) - m - a 2 



For zero doping (SfJ, = 0), for SC, A 
we have 



(17) 



An and 



ImXo 



Nn 



1 



(18) 



Jmxo (z) = No- 



120) 



and 



Im X l ++ {z) = N ^z 



(21) 



So the main result is the same as with no doping - 
Im Xo ± ( z ) i s singular at z = 0, the singularity in Im 
Xo ± (0) gives rise to the singularity in Re Xo ++ (0) an d 
to resonance at a smaller frequency. On the other hand, 
Im Xo ++ ( z ) i s non-analytic at z = 0, but not singular, 
hence Re Xo ++ (0) i s a ^ so n °t singular, and Im Xo ++ ( z ) 
remains non-singular below the bottom of the particle- 
hole continuum. We further confirm these statements by 
performing the numerical calculations for our four-band 
model. We show the results in FigfTOl 
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